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Abstract. The dynamics of a close binary system of globular clusters is considered. 
It is shown that the star transfer process from one of the components to the other 
should lead to the decrease of dimension of the first cluster with simultaneous 
increase of the dimension of the second cluster. The rates of variations of clusters' 
sizes due to the mutual tidal interaction between the components, as well as of 
the orbital rotation of the binary system on the evolution of clusters' sizes, are 
estimated. The results may be attributed to the double nuclei observed in certain 
galaxies. 

1 The statement of the Problem 

The contemporary views on the evolution of globular clusters include the 
following effects: 

i. Evaporation of stars. 

ii. The radius (tidal) of a cluster is determined by the condition of com- 
pensation of the cluster's gravity and the field of the Galaxy. 

iii. Core collapse. 

iv. Core collapse may be halted by one or several hard binaries. 

v. Tidal interaction or shocks of the Galaxy can accelerate the the core 
collapse. 



The time scales of these processes for typical globular clusters are of the 
order of billion years. 

Consider now a binary system both components of which are globular 
clusters, i.e. when we have a binary globular cluster, particulary a cose 
binary globular cluster. Our problem here is to reveal heuristically the prin- 
cipal effects which arise when two ordinary globular clusters can be compo- 
nents of a binary system. 

The problem of binary star clusters can have direct relation to the galax- 
ies with double nuclei, if assuming that its both components are supermas- 
sive star clusters. 



Consider a binary system with globular clusters A and B of masses Ma and 
Mb, mutual distance a and orbital period P by relation 



where the summary mass M = Ma + Mb is in units of solar mass, P in 
days and a in pc. For a binary of globular clusters, for example, with 
M = 1O 6 M , and intercomponent distance a = 300 pc, from (1) we have 
for the orbital period: P = 5 ■ 10 8 years, i.e. comparable with the dynamical 
time of the galaxy. During that time scale such a binary system can perform 
two-three tens of revolutions around the mass center of the galaxy. 

However, by reasons which we discuss later, we shall deal with masses 
of the order of 10 s — 10 9 M Q , i.e. higher as compared with the mass of 
ordinary globular clusters in our Galaxy. Adopting therefore, M = 10 8 M 
and a = 300 pc we find from (1): P = 5 • 10 7 years, i.e. shorter as compared 
with the dynamical time scale of the galaxy. 

The actual sizes of globular clusters in the Galaxy are essentially deter- 
mined by the tidal field of the Galaxy. For an isolated globular cluster the 
tidal radius r g , is determined by the gravitation action of the Galaxy of a 
mass M g and radius R g ; when the cluster is located in the environments 
of the Galaxy the tidal radius is given by so called Hoerner's formula (von 
Hoerner, 1957): 



2 Tidal Radii of Clusters 



a = 9.49 • 10~ 8 P 2 I Z M X I Z , 



(1) 





As we see, the tidal radius r g of a globular cluster depends on the masses 
both of cluster Ma and the Galaxy M g . At Ma = 10 7 M Q , for example, for 
M g = 10 11 M© and R g = 15 000 pc one has r g = 550 pc. 

The situation is quite different when we deal with a close binary globular 
cluster. In this case, the tidal radius r t for one of the components - the 
cluster A, is determined by the gravitation action of the second cluster B: 

GM B GM b GM a 



(a — r t ) 2 a 2 r. 



2 ~ .2 (3) 



or 



1 = 1, (4) 



where x and q denote 



(1 — x) 2 x 

x = r i- q= it- (5) 

Eq.(4) is the correct relationship for the determination of x, i.e. of rt- In 
particular case, when x « 1, which corresponds "Galaxy - Globular clus- 
ter" combination, we shall have from (4): x = (g/2) 1 / 3 or r t = R g (q/2) 1 / 3 , 
i.e. the Hoerner's formula (2). 

As follows from this relationship the tidal radius r t for the component A 
depends on the ratio of masses of both components, but not on the absolute 
magnitudes of their masses. 

In Table I, the numerical values of tidal radius r t , obtained with the help 
of Eq.(4), are presented for three values q = 1, 2, 3 and the intercomponent 
distance a = 500 pc. As we see, in this case the tidal radii are 1.5-2 times 
smaller than those determined by the gravitational action of the Galaxy. 
Particularly, for equal masses of both cluster-components, i.e. when q = 1, 
we have for the tidal radius r t = 265 pc essentially smaller than we have at 
"Galaxy - Globular cluster" combination, i.e. r g = 550 pc. 

Thus, we arrive at the first remarkable, and important conclusion: the 
tidal radius of a globular cluster being a member of a binary cluster system, 
is determined by the gravitation attraction of the second cluster and not of 
the Galaxy. Such clusters, therefore, are expected to be more compact than 
typical isolated globular clusters. In this case the real sizes of globular clus- 
ters are determined not by tidal interaction of Galaxy mass but by mutual 
tidal interactions of globular clusters. 



Table 



I 



Tidal radius r t (A) for the component A of a binary globular cluster 
determined by the gravitation action of the component B at three 
values of mass ratio q and intercomponent distance a = 500 pc 



x 



n(A), pc 



1 



0.53 



265 



2 



0.60 



300 



3 



0.65 



325 



3 Roundchrom of a Binary Globular Cluster? 

Recently the concept of stellar roundchroms for close binary stellar systems 
has been proposed to explain the observed anomalously high luminosities 
of the ultraviolet doublet 2800 Mgll of RS CVn type systems (Gurzadyan, 
1996, 1997 ab). Physically the roundchrom is a common chromosphere 
enveloping both components of the binary system having no geometrical 
contact of their photospheres. The peculiarity of that concept is in the pos- 
sibility to provide a large emission volume around the secondary component 
of the binary system. Geometrical and physical parameters of such config- 
urations have been revealed for over 60 RS CVn type close binary stellar 
systems. 

That idea had inspired us to look for an associated phenomenon for a 
binary globular cluster. As an illustration, in Fig. 1 the configuration of a 
cluster roundchrom, i.e. an eight-shaped zero-velocity equipotential curve 
enveloping both components of the binary cluster is shown for the masses 
of components Ma = 4 • 10 9 M Q , Mb = 2 • 10 9 M Q , i.e. for the mass ratio 
li = 0.33, intercomponent distance a = 800 pc and Jacobi constant C = 3.94, 
so that a narrow corridor can be formed in the intermediate Lagrangian point 
L\ between components. The dashed circles denote the Roche lobes with 
radii i?L determined according to the relationship (Eggleton, 1983) 



0.49q 2 / 3 



(6) 



0.6g 2 /3 + l n (i + g i/3) 



yielding R L {A) = 350 pc and R L (B) = 175 pc. 



In the case of stars, owing to the existence of strong outer boundaries - 
photospheres, one can compare the star's size with its Roche lobe. In ma- 
jority of cases the sizes of main components in close binary systems coincide 
with their Roche lobes. Such a coincidence, again as a rule, is not in the 
case of secondary components; thus an extensive emission volume is formed 
between the outer boundaries of stars and their Roche lobes. 

Fig. 1 

In the case of globular clusters the situation is obviously quite differ- 
ent since one deals with N-body system and hence no clear boundary of 
the system exists. Nevertheless, the chaotic properties of the stellar orbits 
(see Gurzadyan, Pfenniger, 1994) ensure that basic effects can arise also for 
binary globular clusters. Therefore, it appears possible to develop the prin- 
cipal idea of a roundchrom for the case of binary globular cluster supposing 
for simplicity that the linear sizes of the component-clusters are proportional 
to and N-&, and the outer boundary of the first component is coincid- 
ing with its Roche lobe, and is somewhat smaller in the case of the second 
component. 

Then Figure 1 should be interpreted in the following manner. If the 
outer boundary of the component A is not far from the point L\, then one 
should expect a transfer or flow of stars from A to B. The surrounding 
volume, i.e. the "roundchrom" of the secondary will be occupied mainly by 
these stars, resulting a monotonous decrease of the mass of the primary A 
with corresponding increase of the mass of the secondary B. As a result, the 
balance between the kinetic and potential energies cannot remain frozen for 
each of the clusters, along with the secular variation of their linear sizes. 

4 The Expansion and Contraction of Clusters at 
Star Transfer Process 

We now aim to determine the rate of the expansion of the second component- 
cluster (B) of the binary system, and the rate of contraction of the first one 
(A). More definitely, we hope to deduce the law of the growth of the radius 
-Rb(^) °f cluster B and the law of the decreasing of the radius -Ra(^) by time, 
i.e. during the star-transfer process from the first component to the second 
one. 

It is convenient to start from the following basic dynamical principles. 



When the moment of inertia of the system can be considered as constant 
within some time interval, i.e. its variation is slow enough, one readily comes 
to the virial theorem 

U(t) + 2T(t) = 0, (7) 

Note, that this relation is obtained not at the infinite time limit, but via 
averaging within finite time interval. 

Considering the slow stationary star transfer process with conservation 
of the moment of inertia (if the clusters, for example, are far from the 
last phases of core collapse), we can assume that the energy conservation 
condition 

T + U = T(t) + U(t), (8) 

is fulfilled for each component cluster separately. Similarly, it is well known 
that for an isolated cluster though the evaporated stars carry away some en- 
ergy, the total energy of the cluster still can be considered constant within 
some long enough time interval while estimating the dynamical characteris- 
tics of the cluster. 

Then, we have the well known relationship between radius R, velocity 
dispersion v Q and the mass M = mN of the cluster 

o G mN , . 

where N is the number of stars, m is the mass of a star. 

Now we can obtain the variations of the sizes for each of these compo- 
nents of binary cluster A and B during of star transfer process through the 
Lagrangian point L\ from the component A to B. 

First consider the behavior of the component B. For the initial values of 
kinetic T Q and potential energy U a we have 

To = N B ^° (10) 

where R Q is the initial radius of the cluster B. 

As a result of star transfer process the cluster B gains an additional 
kinetic energy 

AT=nt^, (12) 



where is the velocity of stars in the transit point L\, and n is the flux of 
stars with a constant rate n. Thus the total kinetic energy of cluster B at 
the moment t will be 



v 2 

Tit) = T + AT = N B ^+nt^. (13) 

Correspondingly we shall have for the potential energy U (t) of the cluster 
B at the moment t 

m = -1 m2{N X nt)2 > < M > 

where Rt is the modified radius of the cluster B. 

The condition which we will use is that the virial equilibrium remains 
valid during this process, i.e. for moment t. This condition has to be safely 
fulfilled during the stationary star transfer process due to violent relaxation 
effects occuring within dynamical (crossing) time scale. Then, substituting 
(13) and (14) into the virial expression (7) and having in view (9), we obtain 
for the law of the variation of the radius Rt{B) of cluster B by time 

MB) (1 + nt/A B ) 2 
Bo l + (v*/v ) 2 nt/N B - { ' 

This same equation can be obtained using the condition of conservation of 
total energy of each of these clusters during the star transfer process. 

Strictly speaking, v* must be slightly larger than vq. However, adopting 
as a first approximation v* = vq, we obtain 

where the stellar flux rate n is considered constant according to the assump- 
tion on the stationarity of the star transfer process. 

In contrast with cluster B, the cluster A contracts during the stellar 
transfer process, i.e. its radius should be decreased. Similar relationship 
can be derived for this C8.SC dbS well, i.e. for the law of decrease of the radius 
i?t(A) of cluster A 

^)=l-JL t . (17) 
Ro N A { ' 

However, the star transfer process from the cluster A to the cluster B 
cannot proceed stationary too long: with the decrease of the diameter of the 
cluster A, its outer boundary will become more and more farther from the 



point L\. Hence, the initially assumed constant star transfer rate has to be 
substituted by a decreasing function 



n t = n e-^, (18) 

where 7 in yr~l characterizes the rate of the variations of the star transfer 
process from A to B. 

In this case the total number of stars n{t) transferred from A to B at 
the moment t, after starting the transfer process (t = 0), will be 

n(t) = — (1 - e-*) . (19) 
7 

As a result, we shall have for the evolution, i.e. for the growth of the 
radius of the cluster B at the moment t the following relationship, again 
assuming = Vq: 

^± = 1 + J£- (1 - e -*) (20) 

HO 7iVB 

and for the decrease of the radius of A: 

^^-wJ 1 -'^ (21) 

Above, at the deducation of Eqs. (17) and (21), two extremely case for 
the rate of star transfer process are considered, namely, with a constant rate 
and with an exponetially decreasing rate. The real situation, we believe, 
should be expected between these two limiting cases. 

In Fig.2 we give the curves of time variations for both radii, -Ra(£) and 
.Rb(*)> f° r the second case, i.e. at exponetially decreasing star transfer rate 
given by expression (18), for two sets of initial parameters, namely, A^a = 
5-10 9 , Ar = 2-10 9 , 7 = 2-10" 9 yr' 1 and n = 1, n = 2 and n = 3 str.yr" 1 
in the first case (at left), and N A = 5 • 10 s , N B = 2 • 10 8 , 7 = 1(T 8 yr" 1 , and 
no = 0.5, no = 1 and no = 2 str.yr -1 in second (at right). 

Fig. 2 

The main difference between both Figures is the time scale, 10 9 years 
in first case and 10 s years, in the second. While the character itself of the 
evolution of sizes of binary cluster components, the contraction in one case 
and expansion in the other, is independent on rate of the star transfer. Note, 
that the growth of the radius of the cluster B takes place faster than the 



decrease of the radius of A. At the end of evolution, formally at t — ► oo, the 
radius of B is doubled while of A is decreased only on 1/5-th of its initial 
radius. 

Also, the volume of roundchrom around component B is decreased mono- 
tonocally during the evolution of the system. The volume of the roundchrom 
around A tests an increasing as well however not so strongly. 

It is interesting to note that during the essential evolution of sizes of 
cluster A, i.e. during of 2.10 9 years, the distance between its outer boundary 
and the point l\ is remained practically unchanged. This important property 
should be explained also by the slowly motion of the Lagrangian point L\ 
to the direction of the component A synchronously with the decreasing of 
its mass. 

Despite of essential differeces in the time scale, in initial total numbers 
of stars in clusters as well in the rate of star transfer process, the main 
character of curves in Fig.2 is nearly the same. This means that we cannot 
say anything for example about the age of binary system or the phase of 
evolution or star transfer process of the pair judging only by direct image 
of that or another interacting binary system in the centre of some galaxies. 

As an illustration, an evolution sequence of the sizes of a binary cluster 
is shown in Fig.3. The initial data at t = are as follows: Na = 5 • 10 9 , 

Fig. 3 

= 2 • 10 9 with the same mass (solar) of stars in both cases, i.e. 
m = Mq, the ratio of initial radii of both components Ra/Rb = 2.5. The 
configurations of roundchroms correspond to star transfer rate 2 star.yr" 1 
and 7 = 2- 1CT 9 yr -1 . The second configuration from above corresponds to 
t = 0.5 ■ 10 9 yr, the third one - to t = 2 • 10 9 yr. 

In Fig.3, the roundchrom-cluster concentric configurations allow us to 
predict the possible existence of the binary globular clusters, particularly in 
centres of some galaxies, one of components of which in two-envelope form, 
i.e. with an outer low density ring around the central dense spheric disk. 

Returning to the star transfer problem, one should outline the following 
aspects as well: 

a. One can always find wandering stars in the vicinity of the corridor in 
Li, ready for the passage through this corridor. 

b. The gravitation attraction of cluster B, should support the accumu- 
lation of individual stars-members of the cluster A in the vicinity of the 



corridor in the point L±, 



5 The Tidal Effect 

Independently to the star transfer process from one component to another, 
each star, say, in the cluster A, individually undergoes gravitation attraction 
by the cluster B. As a result, the star may acquire an additional velocity. If 
so, the macrostructure of energy balance of cluster A must be modified due 
to the variation of its linear size. 

To evaluate quantitatively these variations, one has to estimate the 
amount of the additional kinetic energy gained due to the gravitation at- 
traction of the total mass Mr of B located in its center on the distance a 
from the center of A as is shown in Fig.4. 

Fig. 4 

Assuming stars of equal mass m homogeneously distributed within 
a spherical cluster A with initial radius i?A and concentration 

we can write down the equation of the motion of a star at the distance r 
from Mb 

dAv mM B 

m ^r = - G —> (23) 

which gives for the additional velocity Av of the star 

Mr 

Av = G^t. (24) 

In view of the gravitational action of the cluster B on individual stars 
in cluster A, i.e. due to the tidal effect, we can write for the elementary 
kinetic energy dT^ of a mass included into the element of sphere of radius 
r and thickness dr 

Av 2 t 2 
dT tid = dm = 7rn mG 2 M| (1 - cos ip) — dr, (25) 

_■ V 

where Mb = m Nb and 



After integration of (25) over limits a — R t and a + Rt, we obtain for 
the additional kinetic energy acquired by the cluster A under the action of 
cluster B 



Ttid(t) " 4 ° m "RT l^ 3 ^ " ^ ^1 ' ' ( } 

where u = a/Rt- 

Then the energy conservation law should be written in the form 

T + U = T(t) + U(t)+T tid (t), (28) 

where To and Uq are the same as above. For T(t) and U(t) we have 

W-f=^. (30) 

Substituting (27), (29) and (30) into (28), and in view of (10) and (11), 
we obtain finally the radius Rt of cluster A at the moment of time t: 

W * fl-*)- Grr, 1« ^ * 2 (31) 



R V ^0/ i? t 4 V^ 2 -l 2u u-1 

with u = a/Rt, and i?o is the initial radius of cluster A. 

The relationship (31) defines the law of the variation of cluster's radius 
R t with time t. As we see, the tidal effect of the cluster B leads to the 
increase of the radius Rt of cluster A. 

In Fig.5 the curves of the growth of cluster's relative radius, Rt/Ro, with 
t calculated with the help of (31), are drawn for three numerical values of 
N A : 4 • 10 8 , 3 • 10 8 and 2 • 10 8 with N B = 10 s and m = M Q , cluster radius 
Rq = 100 pc and distance between the centers of clusters a = 300 pc. 

Fig. 5 

Note, that the growth of the radius Rt/Ro of cluster A due to the tidal 
effect of the cluster B, occurs rather rapidly - during a time scale of the 
order of 10 7 yrs. This result seems to be in disagreement with the dynamical 
nature of the globular clusters. However, as we will show in the next section, 



the orbital rotation effect can essentially suppress the rapid growth of the 
radius of cluster A. 

The peculiarity of the problem considered above, namely, the role of the 
tidal process in binary globular clusters, is that one deals not with atoms 
in the stellar atmospheres of interacting close binaries but with stars. The 
main difference is in the absence of collisions between stars, i.e. the absence 
of energy exchange process between stars (particles). 

The next step may be, perhaps, the examination of both processes, the 
tidal and star-transfer, simultaneously. However, as shows the preliminary 
analysis, the results and basic conclusions will remain unchanged. 



Each star in both clusters will undergo the action of centrifugal force pro- 
voked by the orbital motion the components. This force may compensate in 
certain degree the gravitational attraction of the other cluster. This factor, 
as appears, may change radically the above conclusions. 

The centrifugal force depends on the tangential velocity V% and the dis- 
tance to the center of rotation of the binary system. Gravitational attraction 
provoked by cluster B can be compensated completely by centrifugal force 
only if the tangential velocity Vt of a star in the center of cluster A satisfies 
the condition 



where Ma and Mb are the masses of the clusters, as before. 

On the other hand, for the orbital velocity V or of a star we have ( P is 
the orbital period of the binary system) 



or using the known expression for P, 



6 The Role of Orbital Rotation 




(32) 




(34) 



From (34) and (32) we find for the ratio of both velocities 



^ = 1(1 + M a /Mb). 



(35) 



In particular case when Ma = Mb, we achieve to an important result 



V t = V 



i.e. only at equal masses of both clusters we shall have an equilibrium of 
both types of forces - centrifugal and gravitational; in this case the tidal 
acceleration of a star towards the cluster A due to the gravitation attrac- 
tion of cluster B will be compensated by the centrifugal force of an orbital 
rotation of the system. The further evolution of the sizes of clusters will be 
determined by star transfer process from one component to another. 

Note, that depending on the numerical value of the ratio Ma/Mb we 
shall have different relations between Vt and V or . So, if Ma/Mb > 1 we 
have 



In the first case, when Vt < V or , the centrifugal force will be compensated 
by the gravitation attraction not completely, hence, in this case the tidal 
acceleration due to the component B may have some role. Accordingly, 
some increase of total kinetic energy or decrease of potential energy in the 
system A and, correspondingly, a decrease of the linear sizes of cluster A 
may be expected, although on less degree as compared with the limiting 
case given by relationship (21) or curves in Fig.2. 

In the second case when Vt > V OT , the centrifugal force will dominate 
over the gravitational attraction, and the final result will be the same as 
in the first case, i.e. we shall have some increase of kinetic energy or some 
decrease in linear sizes of the cluster. 

The condition Vt = V or can hardly be fulfilled for real cluster binaries, 
therefore the general conclusion may be formulated in the following manner: 
the decrease of linear sizes of the main component of a binary cluster at star 
transfer process is inevitable although at a lower rate as compared with the 
rate given by (31). 

As to the rate of variation of cluster linear size, it depends in which 
degree the additional velocity Vt is larger or smaller than the "thermal" 
(mean) velocity of stars Vq given by relationship 



V t <V 



and vice versa, when Ma /Mb < 1, we shall have 



V t >V c 



or 




1 GM A 



(36) 



Prom (32) and (36) we find 



V 




1/2 



(37) 



For example, at M A /M B = 2.5 and a/R A = 3 we have: V t /V = 0.28, 
i.e. not too small in order to be ignored and not too large in order to be 
compared with the results presented in Fig.4. 

Thus, the star-transfer process from cluster A to cluster B leads to the 
decrease of the size of cluster A. The tidal effect of the cluster B on cluster 
A leads to the growth of the sizes of A. These are the consequences of 
the tidal effect in its "pure" form. However, the rotational effect due to the 
orbital motion of the clusters A and B, essentially reduces the rate of this 
process. As a result, the behavior of the size of component A has to be 
rather complex. 

Thus the star-transfer process from cluster A will definitely be accompa- 
nied by the decrease of the cluster's size, however with a rate essentially de- 
pended on various physical, dynamical and kinematical parameters of both 
clusters. 

The latter two effects, the tidal and orbital rotation, lead to changes in 
dynamical state of globular cluster during a time scale of the order of tens 
of millions years, which is much shorter as compared with the time scale of 
globular cluster's evolutionary processes mentioned in Section 1. Thus the 
clusters in the binary systems must have different evolutionary time scales 
than the isolated globular clusters. 

7 On the Binary Nuclei of Galaxies 

During the last decade the discoveries of galaxies with double nuclei become 
more and more common. Due to their form and structure the nuclei are 
usually interpreted as binary black holes. Nevertheless the the concept that 
nuclei are compact stellar systems cannot be excluded, i.e. a supermassive 
globular clusters - of a mass of components of the order of 10 9 M & and more. 
An example of such a double nuclei is shown in Fig.6; this Mrk 273 image 
is taken by Knapen et all. (1997) via Keck telescope. This image resembles 
an eight-shaped binary configuration with different sizes of components. At 
the distance of this galaxy 160 million pc, the intercomponent distance is 
estimated as 800 pc, with diameters of components 450 pc and 700 pc. 



Table 



II 



Total number of stars JVa and Nb in the components A and B of binary 
nuclei of six galaxies with known absolute magnitudes Ma and Mb of 



components 


Mrk, No 


M A 


M B 


N A 


N B 


266 


- 17.8 


- 17.5 


2.1 • 10 9 


1.6- 10 9 


273 


- 18.4 


- 17.7 


3.6 


2.0 


463 


- 19.5 


- 19.3 


10 


8.3 


673 


- 19.6 


- 19.6 


10 


10 


739 


- 19.1 


- 18.3 


6.9 


3.3 


789 


- 19.5 


- 17.5 


10 


1.6 



Fig. 6 

Note that in Fig. 6 the space between the components reveals certain 
brightness. According to our concept the point L\ located there must pro- 
vide a corridor for star flow from one component to another, so that the 
inter component space must acquire definite brightness, which can even be 
roughly estimated. At the flux 1 star per year with stellar velocity 100 
km s _1 we shall have 10 7 stars during 10 6 yrs on 10 pc linear length on this 
area, thus forming a noticeable background brightness, as maybe indicates 
Mrk 273. 

At present the number of known galaxies with binary nuclei is over one 
hundred. For some of them, including Mrk 273, the absolute magnitudes of 
the components Ma and Mb have been obtained (Khachikian, 1998). Then, 
assuming for these nuclei completely stellar composition of solar mass and 
solar luminosity, we can evaluate the total number of stars A^a and Ab in 
the components; the results are presented in Table II. 

In the case of Mrk 273, the total number of stars in the components is 
found 3.6 • 10 9 and 2.0 • 10 9 , i.e. almost of the same order as those used 
above when drawing the eight-shape curve in Fig.l. Moreover, the derived 
ratio Na/Nb ~ 2 for this galaxy corresponds, at least qualitatively, to the 
photographic images of components in Fig. 6. Note, the existence of nuclei 
with equal number of stars in the components (Mrk 673, ratio 1:1) as well 



with the ratio up to 1:5 (Mrk 789). It is remarkable that the total number of 
stars in the binary nuclei at least for these six galaxies is of the same order 
- 10 9 -10 10 . Even in view of observational selection, this fact cannot be 
totally ignored concerning the nature and composition of the binary nuclei 
of galaxies. 

Fig. 7 

Another example we have in the case of IC 4553 (Arp 220) with a re- 
markable binary nuclei structure discovered at radio continuum 4.83 GHz, 
see Fig. 7 (Baan, Hischick, 1995). The projected separation between the 
components is 330 pc, the reconstructed separation 466 pc. At orbital period 
of 7 • 10 6 yr for the two nuclei with equal masses, the estimated dynamical 
mass for both components is 10 10 Mq. It is important to note, that in this 
case no clear hot spots are found in the centers of components within the 
line emission structure. The spectral properties confirm the starburst domi- 
nated nature of the observed emission, and hence, the stellar content of both 
nuclei. 

The results obtained in present article may find an interesting appli- 
cation, particularly, for binary globular clusters in the Large Magellanic 
Clusters. 

8 Conclusions 

Thus, we formulated and qualitatively considered the binary system with 
globular clusters as components reveals the importance at least of three as- 
pects on their dynamics: 

- Mutual tidal interaction of both components. 

- The role of the orbital rotation of the binary system. 

- The star transfer process from one component of the system to the 
other. 

The preliminary conclusions can be summarized as follows: 

• Globular cluster as a component of a binary system gains new dynam- 
ical characteristics and properties; 

• Dynamical evolution of a globular cluster as a member of binary system 
occurs more rapidly as compared with an isolated cluster; 

• The tidal radius of a globular cluster in a binary system is smaller as 
compared with the tidal radius determined by the Galactic field; 



• The star transfer process from one component to another in a binary 
cluster system leads to the reduction of linear size of the first component 
and the increase of the size of second one; 

• In binary globular clusters the tidal interaction leads to the increase 
of their kinetic energy and, hence, to the expansion of clusters. 

• Orbital rotation of a binary cluster system plays a regularing role on 
the variations of the sizes of the clusters. 

Extensive numerical study of this problem can be of remarkable interest. 
I thank V.G.Gurzadyan for valuable discussions. 
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Figure Captions 

Fig.l. Roundchrom structure for a close binary system with star clusters as com- 
ponents with mass ratio \i = 0.33, masses of components Ma = 4 • 10 9 M Q , 
Mb = 2 • 10 9 M Q and intercomponent distance a = 800 pc, -Rl(A) and -Rl(B) 
are the radii of Roche lobes. All linear sizes are in parsecs. 

Fig. 2 The time-dependent curves of the variations of linear radii Rt/Ro of both 
components, A and B, of binary system with globular clusters as components. 
For the initial moment, t — 0, it is assumed Rt/Ro = 1- The curves are 
calculated for three values of star-transfer rate n D from 0.5 up to 3 str.yr -1 



and for two values of gamma, 2.1CP 9 yr -1 and 10~ 8 yr -1 and the time 
scale 10 9 yr and 10 8 yr. During the star-transfer process from component 
A to component B, component A undergoes a decrease of its radius, i.e. 
Rt/Ro < 1 (lower half) and the component B, the contrary, an increase, 
Rt/Ro > 1 (upper half). 

Fig. 3 Three stages of the evolution of globular clusters' diameter as a function of 
time t, i.e. during star-transfer process from component A to component B 
through the corridor in Lagrangian point L\. All three configurations are 
calculated for one and the same star-flow rate n = 2 str.yr -1 and gamma = 
2.10~ 9 yr -1 . Initially, at t = 0, the component B has a powerfull roundchrom 
(pointed area) and component A is practically without a roundchrom. 

Fig. 4 To the calculations of the tidal effect of the cluster B of a mass Mb on the 
stars in cluster A. 

Fig. 5 The estimated time-depended variations of radius, R t /Ro, of component A 
due to the tidal effect for three values of total number of stars Na in the 
cluster A, at its radius R = 100 pc and intercomponent distance a = 300 
pc. 

Fig. 6 Double nuclei of Mrk 273 (Knapen et al 1997). Components of nuclei are of 
different masses and different sizes. The separation between components is 
730 pc. Note the background in the corridor between components. 

Fig.7 Double nuclei of IC 4553 in 4.83 GHz (Baan and Haschick 1995). The 
separation between nuclei is 330 pc. 
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